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Abstract 

The bounded orbital motion of a massive spinless test particle in the 
background of a Kerr Brans-Dicke geometry is analysed in terms of 
worldlines that are auto-parallels of different metric compatible space- 
time connections. In one case the connection is that of Levi-Civita 
with zero-torsion. In the second case the connection has torsion de- 
termined by the gradient of the Brans-Dicke background scalar field. 
The calculations permit one in principle to discriminate between these 
possibilities. 
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1 Introduction 



In general relativity it is commonly assumed that idealised massive spinless 
test particles have spacetime histories that coincide with time-like geodesies 
associated with the spacetime metric when acted on by gravitation alone. 
Such assumptions embody such notions as the "equivalence principle", the 
"equality of inertial and gravitational mass" and the tenets of "special rela- 
tivity" . Such a geometrical description of gravitational interactions is formu- 
lated in terms of Einstein's torsion-free metric-compatible (affine) connection. 
The particle histories then have self-parallel 4- velocities (tangent vectors) and 
may be termed "Levi-Civita auto-parallels" . If a massive test particle moves 
in a bounded orbit in the geometry described by the exterior Schwarzshild 
metric one may use such hypotheses to calculate its orbital perihelion shift 
per revolution. Despite competing perturbations, such a perihelion shift of 
the planet Mercury in the gravitational field of the Sun is regarded as one 
of the classical tests of Einstein's theory of gravitation [1] , [2] . Further ev- 
idence for these hypotheses is sought from observations of the pulse rate of 
the binary pulsar PSR 1913+16 that appears to be speeding up due to grav- 
itational radiation [3]. A number of efforts have been made [4], [5], [6] to 
prove the geodesic hypothesis for test particles from a field theory approach 
but none are entirely convincing. However, (for spinless test particles) this 
hypothesis has almost become elevated to one of the natural laws of physics 
[7]. 

In 1961 Brans and Dicke [8] suggested a modification of Einsteinian grav- 
itation by introducing a single additional scalar field with particular gravita- 
tional couphngs to matter via the spacetime metric. They suggested that the 
experimental detection of such a new scalar component to gravitation might 
have been overlooked in traditional tests of gravitational theories. Their the- 
ory is arguably the simplest modification to Einstein's original description. 
Since then a number of astrophysical observations have indicated that grav- 
itational scalar fields may indeed have relevance to the dynamics of matter. 
Furthermore on the theoretical side, modern "low energy effective string the- 
ories" are replete with scalar fields and most unified theories of the strong 
and clcctroweak interactions predict such fields with astrophysical implica- 
tions [9]. 

Although the gravitational field equations were modified by Brans and 
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Dicke [8], they still assumed that the motion of test particles under the in- 
fluence of gravitation was described by a Levi-Civita auto-parallel associated 
with the metric derived from the Brans-Dicke field equations. Dirac later 
suggested [10] that it was more natural to generate the motion of a test par- 
ticle from a Weyl invariant action principle and that such a motion in general 
differed from a Brans-Dicke Lcvi-Civita auto-parallcl. Although Dirac was 
concerned with the identification of electromagnetisni with aspects of Weyl 
geometry, even for neutral test particles it turns out that such test particles 
would follow auto-parallels of a connection with torsion. In Ref. [11] we have 
shown that the theory of Brans and Dicke can be reformulated as a field the- 
ory on a spacetime with dynamic torsion T determined by the gradient of 
the Brans-Dicke scalar field 

r = e«® — ®e"®X„ (1) 

in any coframe {e"} with dual {X^}. Although (in the absence of spinor 
fields) no new physics of the fields can arise from such a reformulation, the 
behaviour of spinless massive test particles in such a geometry with tor- 
sion is less clear cut. Two natural alternatives present themselves. One 
may assert that their histories are either time-like geodesies associated with 
auto-parallels of the Levi-Civita connection (as did Brans and Dicke) or the 
auto-parallels of the non-Riemannian connection with torsion. In [12, 13] 
we have shown that it is possible to theoretically compare these alternatives 
for the history of a mass in orbit about a spherically symmetric source of 
scalar-tensor gravity by regarding it as a spinless test particle. In principle 
such calculations could be confronted with observation in space experiments 
that measure the orbital parameters in binary systems. In practice there are 
many competing perturbations that may contribute significantly. In particu- 
lar a rotating gravitational source will be expected to produce a modification 
whose magnitude will depend on both its mass and angular momentum as 
well as the scalar field. It is of interest therefore to compare the auto-parallel 
orbits derived from the two connections above for a spinless test particle in 
the metric of a spinning source. 

In this note we perform such a calculation for a particle in the geometry 
derived from the exact Kerr Brans-Dicke metric and scalar field found in [14] 
(See also [15]). This solution to the Brans-Dicke equations is asymptotically 
flat and reduces to the Kerr metric solution to Einstein's equations when the 
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scalar field is replaced by a constant. The solution describes a stationary and 
axially symmetric metric and depends on parameters that may be identified 
with the scalar charge, mass and angular momentum of a localised source. 



2 The motion of massive test particles 



The Kerr Brans-Dicke solution found in [14] can be written in Boyer-Lindquist 
coordinates {t, r, 9, (p) as; 

4A 



r-{M+JM^-P)\ SA „ Psin^^,^ 2M/r ^ 

^ ' \r-{M-^M^-P)l ^ P E P ' 



2kA^ 



„2 



/ {r - MY - {M"^ - \ 2 dr 

with the scalar field, 

^''\r-(M- ./M^) ) ^ ' 

where M > I denotes the source mass and I its angular momentum per 
unit mass. The constant A determines the strength of the scalar field and 
k — cu + 3/2, in terms of the Brans-Dicke parameter cu. As usual 

E = + /2 cos^ e 

A = + /2 _ 2Mr (4) 
P = AE + 2Mr(r2 + i2). 

We first examine test particle orbits C for massive spinless particles that 
follow Levi-Civita autoparallels. The equations of motion are 

VcC = 

in terms of the torsion-free Levi-Civita connection V and 4-velocity C nor- 
malised according to 

g(C,C) = -l, (5) 
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(throughout units are adopted such that G = 1 and c = 1). If C : r i— >■ a;'*(r) 
in terms of proper time r these equations yield 



d ( d/Cc^ \ r 1 1 ^ dec dec „ / _ \ 

The metric above has two independent KiUing vectors dt and d^. These 
generate two constants of motion: the particle energy and orbital angular 
momentum. Since the orbits are planar we take ^ = f and set 

E ^m{igtt + (pgt<f): (7) 

L = m{igt^ + (fg^^). (8) 

One can express r in (5) in terms of t and (f and metric components on the 
orbit. Eliminating t and from equations (7) and (8) and substituting into 
(5), the orbit equation for the particle equation may be written: 



dip J {g^^E - guLYgr 



{$-2$oA + 2gt^EL - guL^ - g^^E^} (9) 



where L = ^^'^'^^ ' and E = Now define, 

m. m. ' 



A{r) = 4M^'^-Ar-^ _2MZr (10) 
B(r) = ^ (11) 



C(r) = ^ (12) 



A 



p^{r) = [f + t^Y + 2Ml\ (14) 
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and introduce u = - so that: 

r 

yd'Pj (b{1/u)E-LA{1/u)) G{1/u) 



C{l/u)E^ -PA{1/u)}. (16) 

To analyse this equation we employ the physically motivated approxima- 
tions discussed in [13]. Thus if the radius of a (weak field) Newtonian orbit 
is much larger than the corresponding Schwarzschild radius of the source, 
one may expand this equation around m = up to third order in order to 
compare its solutions with those in a Schwarzschild background: 

(^) -So + Siu + S2u'' + Ssu' (17) 



where the constants: 



So = i^,{E'-l), (18) 



1 ....^.^o 1 



Si = -r^4MlE{E^ - 1) + ^(2M - VM^-l^A), (19) 



S2 = -l + ^{-l{M^-f)A^ + 3f{E^-l) 
2 

+MVM2 - PA + 2k{M^ - l'^){E^ - 1)A^} 



+ ^{8M^EH - 4MIEVM^ - PA] 



+ 11 (^'-1), (20) 



^3 = 2M + J-{-3ZVM2 - PA - -VM2 - P{AM'^ - f)A 

3 

+6Ml^E^ + 2MVM2 - PA + AkE^M{M^ - l'^)A^ 
-{2k + -) - A^} + \{-2MlE{M^ - /2)A' 



2 



6^ V J ' L3 

-AMHEy/M'^ - PA + 
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+8kEMl{E' - 1)(M' - l')A' + 12MEl\E' - 1)} 
+i^{2M^^2;'(2^' - 1) - M^EH'^VM^ - PA} 

+^32fM^E^{E'^ - 1). (21) 

All the terms in S2 except —1 and all the terms in ^'3 give corrections to the 
Newtonian orbital equation. 

By contrast we now compare this equation with the one obtained by 
assuming the worldline is a timelike auto-parallel of a particular connection 
with torsion: 

VfjC = 0. 

Here V denotes the connection with the torsion (1) specified by the scalar 
field in terms of the 2-forms 

T'^ = e« A — 
2$ 

and the 4-velocity C is again normalized with 

It is possible to re-write the worldline equation in terms of the Levi-Civita 
connection V as 

where for any vector field V, V = g(V^, — ) is the metric related 1-form. This 
may be further simplified to 

restricted to C, which in local coordinates gives: 

d ( ^M')dx^^\ ^Morn , dx" dx-^ „„ dy^ 

_ + = -r^.. (22) 

For any Killing vector K with — 0, the expression 

7K = $'/'g(K,C) 
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is constant along the worldline of the particle. As before the Killing vectors 
Kt = dt and K^^ = generate, respectively, two constants of motion E and 
L, corresponding to energy and orbital angular momentum. Thus 

^ = it9tt + <f9t^), (23) 



^^"^{^I i^9t^ + ^9vv) (24) 



in terms of metric functions evaluated on planar orbits. Eliminating t and 
(p from equations (23) and (24) and substituting in (5), one obtains the new 
auto-parallel orbit equation as 



-2 



{$-iA + 2g,^EL - guL^ - g^^E^} (25) 



where E = '^^^"^ ^ and L — ^liial^ Expressed in terms of the variable 
M = -, it becomes 

h- = 7 : . \2 u) + 2B{l u)EL- 

(b{1/u)E-LA{1/u)) G{l/u) 

C{l/u)E'' -fAil/u)}. (26) 
Expanding to third order in u as before, one finds 

(^) -Co + Ciii + C2ii' + C3ii=^ (27) 



in terms of the constants: 



Co = i(£^^-l) (28) 



C. = 2i^ + 4^(i^-l) (29) 

C2 = j^{3P + 2kA\M^-f)}{E^-l) + i^8MHP 

+i^l2M^l^E^{E'' - 1) - 1 (30) 



+ i{[(12M - SkMA^^ + 8kMHA^]E{E^ - 1) + IGM^E^/} 
+i^{24M^Z2^2(2^2 _ ^ J-{32M^Z3^^(^2 _ (g^^ 

The first three terms of C2 and all terms in C3 imply corrections to Newtonian 
orbits. 

We note that both orbit equations have been written in the form: 

(du \ 
~ ^(m) = Lo + LiM + Ls-u^ + Ls-ul (32) 

so their solutions can be analysed in terms of the corresponding constants 
according to the roots of the equation g{u) — 0. Suppose first that all three 
roots are real, i.e. 

AL\L^ + ALqLI - LlLl + 27LlLl - I8L0L1L2L3 < 0. 

Suppose further that the roots are distinct and ordered to satisfy i^i < ^2 < 
1^3. Then 

L2 

^3 

Prom the orbit equation, g{u) > throughout the motion. Thus, g{u) will 
have a local maximum between ui and U2- Hence, for a bounded orbit, 
Ui corresponds to aphelion and U2 corresponds to perihelion. Consider the 
following cases [16]: 

i. If Ui > 0, one obtains bounded orbits of "elliptic "type. This requires 
that both Lq and L2 be negative provided that L3 > 0. Then the particle 
is confined to the interval ui < u < U2. (If ui = U2, one obtains circular 
orbits.) 

ii. If ui — 0, one obtains open orbits of "parabohc "type. This requires 
that Lq = 0. This is possible for orbits associated with both Levi-Civita and 
torsional connections provided E"^ = $0 ^m^. 

iii. If Ui < 0, one obtains open orbits of "hyperbolic" type. This requires 
that E'^ > ^Q^m^ provided that L3 > where L3 = C3 if the orbit is 
associated with an autoparallel of the torsional connection and provided L3 = 
S3 if it is associated with the Levi-Civita connection. 
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3 The analysis of bounded orbits 



We are interested here in bounded orbits [17] in which case the general so- 
lution of (32) can be expressed in terms of Jacobian elliptic functions. By 
introducing the variables 



X = -(pjLs {us - ui), 



U — Ui 



2 ^ y M2 — Hi 

(32) becomes 



(I) -i^-y'W-'^y^ (33) 



with p = /«2_ui_ Its general solution is 



y ^ sn{x + 5) (34) 

where 5 is an arbitrary constant. Hence for both connections yielding orbits 
with perihelia, 



u 



-Ux^ {U2 - Ui) STl? {^ip^L2.{u2, - Ui) + . (35) 



The periodicity of these solutions enables one to calculate a perihelion shift 
per revolution. The increase in </? between successive perihelia is given pre- 
cisely by 

= 2 / , (36) 

■^"1 ^J Ls{u - Ui){u - U2){u - U3) 

With the transformation y = ./ """"^ . this becomes 

AK 

A^=^=_ (37) 

where 

K 



^(l-y2)(l_p2^2)- 
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Depending on circumstances one may be able to approximate this inte- 
gral. This is possible if one is interested in non-relativistic bounded orbits 
in which the dimensionless quantity Mu remains small compared with unity 
and the orbital speed is small compared with the speed of light. This would, 
for example, arise for the motion of the planet Mercury regarded as a test 
particle in orbit about the Sun as a source. Even at the Sun's surface, where 
Rq = 7 X lO^m, Mq = 1.477 x lO^m , Mu = 2.11 x 10"^. At aphelion 
and perihelion, both —j^Ui and —j^U2 are small quantities [18]. Thus in 
the following we approximate — ff'^a — 1- This means can be considered 
small so: 

^ - ^TT (l + 

and since the ratios — and — are small, we expand 

Us U3 ■) f 

2 {U2 - Ui) {U2 - Ui)Ls 

~ U3 L2 ' 

We further approximate the term 

L 



l-^(2«i + «2)) 



After a little algebra one finds that the increase in (p per revolution becomes 
A^c.^(l-^^iu, + U2)). (38) 



The advance of perihelion (perihelion shift) per revolution would be 

S = Aip- 27T. 

It is conventional to define e, the eccentricity of an elliptic orbit and ro its 
semi-axis major. Then 

{1 + e)ro, r2 = (1 - e)ro 

where 1^1 = :^ corresponds to aphelion distance and r2 = corresponds to 
perihelion distance of an elliptic orbit. The perihelion shift may be expressed 
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in terms of these orbit parameters for the Hmiting Newtonian Kepler elhpse. 
Once a set of Kepler orbit parameters have been ascertained then these for- 
mulae permit one to match them to a relativistic orbit in terms oi M, I, cu 
and A and the constants of motion [13] . 

Thus in summary, we have examined the bounded orbits of test particles 
in the presence of a particular solution to the Brans-Dickc theory of gravita- 
tion. Even with the same constants of motion and the same limiting Kepler 
orbits such orbits will in general differ. Those with histories determined by 
the Levi-Civita connection and satisfying the approximations above, have a 
perihelion shift per revolution 



By contrast if the history is determined by the connection with torsion, the 
perihelion shift per revolution is different, being given by 

^/\C2\\ 2(l-e2)ro|C2|y 

Finally we note that, when A = (i.e. the scalar field $ is constant) both 
orbit equations describe geodesic motion given by the Levi-Civita connection 
in a background Kerr geometry. If one further sets I — 0, they both describe 
geodesic motion in a background Schwarzschild geometry. In this case the 
constants reduce to C3 = S3 = 2M and C2 — S2 — —1 and the perihelion 
shift reduces to the classical value [18] 

GttM 



4 Conclusion 



An analysis of bounded orbital motion of a massive spinless test particle in 
the background of a Kerr Brans-Dicke geometry has been given. Such orbits 
have been discussed in terms of worldlines that are auto-parallels of different 
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metric compatible space-time connections. In one case the connection is 
torsion-free while in the other, the connection has a torsion that depends on 
the (spacetime) gradient of the scalar field. Such a connection differs from 
that of Levi-Civita by terms that contribute to the so called "improved stress- 
energy tensor" that enters into the original formulation of the Brans-Dicke 
theory. The variational derivation of this theory in terms of independent 
metric and connection variations naturally gives rise to a connection with 
such torsion. 

It is notoriously difficult to discriminate between gravitational theories 
by observing planetary orbits in the solar system since many perturbations 
need to be taken into account. In the above calculations one has some infor- 
mation about the value of the constant u) but little guidance as to the value 
of the constants $o s-nd A. Furthermore there is no unique Kerr Brans-Dicke 
geometry that one may, with confidence, ascribe to any particular spinning 
source. Despite these shortcomings, we feel that there may be astrophysical 
situations where considerations relevant to non-Riemannian geometries be- 
come relevant. The motion of bodies around spinning black holes, neutron 
or magnetic stars are examples. In order to obtain evidence for any grav- 
itational interaction one needs to analyse observational data in terms of a 
particular theoretical framework. We have suggested that the detection of 
scalar-tensor gravitational interactions may benefit from a non-Riemmannian 
description. The calculations presented here offer such a framework in the 
context of the Brans-Dicke geometry with and without torsion. The effects 
of coupling other types of fields such as a Kalb-Ramond 2-form potential 
maybe worth investigating [19]. 
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